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Abstract: Photonic crystal nanobeam cavities are versatile platforms of 
interest for optical communications, optomechanics, optofiuidics, cavity 
QED, etc. In a previous work (TJ, we proposed a deterministic method 
to achieve ultrahigh Q cavities. This follow-up work provides systematic 
analysis and verifications of the deterministic design recipe and further 
extends the discussion to air-mode cavities. We demonstrate designs of 
dielectric-mode and air-mode cavities with Q > 10 9 , as well as cavities with 
both high-Q (> 10 7 ) and high on-resonance transmissions (T > 95%). 
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1. Introduction 

High quality factor (Q), small mode volume (V)[2 | optical cavities provide powerful means for 
modifying the interactions between light and matterO, and have many exciting applications 
including quantum information processing!!], nonlinear opticsQ, optomechanics[6|, optical 



Fig. 1. (a) Schematic of the proposed nanobeam cavity, (b) FDTD simulation of the energy 
density distribution in the middle plane of the nanobeam cavity. 

trapping[7] and optofluidics[8|. Photonic crystal cavities (PhC)ll9l lfT0l have demonstrated nu- 
merous advantages over other cavity geometries due to their wavelength-scale mode volumes 
and over-million g-factors[ 1 1 1 - [ 30 1 . Although small mode volumes of PhC cavities can be eas- 
ily achieved by design, ultrahigh Q factors are typically obtained using extensive parameter 
search and optimization. In a previous workjl], we proposed a deterministic method to design 
an ultrahigh Q PhC nanobeam cavity and verified our designs experimentally. The proposed 
method does not rely on any trial-and-error based parameter search and does not require any 
hole shifting, re-sizing and overall cavity re-scaling. The key design rules we proposed that re- 
sult in ultrahigh Q cavities are (i) zero cavity length (L = 0), (ii) constant length of each mirror 
Cperiod'=a) and (iii) a Gaussian-type of field attenuation profile, provided by linear increase in 
the mirror strength. 

In this follow-up work, we provided numerical proof of the proposed principles, and system- 
atically optimized the design recipe to realize a radiation limited cavity and waveguide coupled 
cavity respectively. Furthermore, we extended the recipe to the design of air-mode cavities, 
whose optical energies are concentrated in the low-index region of the structure. 

Nanobeam cavities have recently emerged as a powerful alternative to the slab-based 2-D 
PhC cavities lfTTl - El . Nanobeams can achieve Qs on par with those found in slab-based ge- 
ometries, but in much smaller footprints, and are the most natural geometries for integration 
with waveguides[23|-[30|. Our deterministically designed cavities have similar structures to 
the mode-gap cavity proposed by Notomi et al. 11231 . and later demonstrated experimentally by 
Kuramochi et al.|l3T|, as well as our own work f32l . We note that the same design principle dis- 
cussed here could be directly applied to realize ultra-high Q cavities based on dielectric stacks 
that are of interest for realization of vertical-cavity surface emitting lasers (VCSELs) and sharp 
filters. Finally, it is important to emphasize that while our method is based on the framework 
of Fourier space analysis|33][34][35], alternative approach, based on phase-matching between 
different mirror segments, could also be used to guide our design, as well as to explain the 
origin of deterministic ultra-high Q-factors in our devices ll36l [371 . 

2. Numerical Verification of the Deterministic Design Approach 

In this section, we use finite-difference time-domain (FDTD) simulations to systematically 
study the design principles proposed in We consider a TE-polarized dielectric-mode cav- 
ity, i.e. Ey,H x ,H z are the major field components, and the energy density is concentrated in 
high index material (silicon in our case). Fig. HI a) shows the schematics of the proposed cavity 
structure H]. It consists of an array of air-holes in decreasing radii, etched into a ridge waveg- 
uide. The hole-to-hole distances are constant ("periodicity"). The structure is symmetric with 
respect to the dashed line in Fig.[TJa). In contrast to the majority of other cavity designs, cur- 
rent structure has no additional cavity length inserted between the two mirrors (L=0), that is 
the distance between the two central holes is the same as that of the rest of the structure (a). 
This design minimizes the cavity loss and the mode volume simultaneously. The cavity loss is 
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Fig. 2. (a) Mirror strengths of each mirror segment under different tapering profiles ob- 
tained from the plane wave expansion method (' 1 ' indicates the mirror segment in the center 
of the cavity), (b) Band diagram of the TE-like mode for / = 0.2 and / = 0. 1 . The green line 
indicates the light line. The circle indicates the target cavity resonant frequency, (c) Mir- 
ror strengths at different filling fractions from the 3D band diagram simulation, (d) Mirror 
strengths as a function of mirror number after quadratic tapering, (e) Radiation-<2 factors 
when nanobeam cavities have different cavity lengths between the two Gaussian mirrors, 
(f) Resonances of the cavities that have different total number of mirror pair segments in 
the Gaussian mirror, and their deviations from the dielectric band-edge of the central mirror 
segment, obtained from FDTD simulation and perturbation theory, (g) H z field distribution 
on the surface right above the cavity, obtained from 3D FDTD simulation. The structure has 
dimension of a = 0.33fim,b = 0.7 fim, the first 20 mirror segments (counted from the cen- 
ter) have / varies from 0.2 to 0.1, followed by 10 additional mirror segments with / = 0.1. 
(h) H : field distribution on the surface right above the cavity, obtained from the analytical 
formula H z = sin(|x)exp(-(jx 2 )exp(-^y 2 ), with a = 0.33/J.m,a = 0.14,<^ = 14. (i) H z 
field distribution along the dashed line in (g)&(h). Length unit in (g)-(i) is fjm. 



composed of the radiation loss into the free space (characterized by Qrad) an d the coupling loss 
to the feeding waveguide (Q W g). Qwg can be increased simply by adding more gratings along 
the waveguide, gi-ad can be increased by minimizing the spatial Fourier harmonics of the cavity 
mode inside the lightcone. The optical energy is concentrated in the dielectric region in the 
middle of the cavity (Fig.[TJb)). In order to achieve a Gaussian-type attenuation, we proposed 
to use a linearly increasing mirror strength along the waveguide) 1|, which was achieved by ta- 
pering the hole radii. To analyze the ideal tapering profile, in this work, we first use plane wave 
expansion method and then verify the results with 3D FDTD simulations. The dielectric profile 
of the structure in the middle plane of the cavity can be expressed as 

'i--^)s(p) (i) 



e(p) £si £s 
with 

s(o)-l 1 \ p - r ^ R 

Tj = j -ax, a is the period, and j = ±1,±2... are integers. R is the radius of the hole. Using 
plain wave expansion method[38 | in the beam direction (i), 

-L = Kb + K x e iGx + K-ie- iGx + ... (2) 
e(x) 

where G = 2k /a. The zeroth (Kq) and first (K"i) order Fourier components can be expressed as 



1 - f 

Ko = f+-^- (3) 

J\ is the first order Bessel function. Filling fraction / = kR 2 jab is the ratio of the area of the 
air-hole to the area of the unit cell. 

The dispersion relation can be obtained by solving the master equation l39l : 

c 2 d 2 E d 2 E 

w ) 



e(x) dx 2 dt 2 

Inside the bandgap, the wavevector (k) for a given frequency (co) is a complex number, 
whose imaginary part denotes the mirror strength (7). For solutions near the band-edge, of 
interest for high-g cavity design|fl], the wavevector k = (1 +ij)n/a, and frequency ft) = 
(1 — 8)y/Konc/a. Substituting this into the master equation, we obtained 8 2 + J 2 = k 2 /4-Kq. 
The cavity resonance asymptotes to the dielectric band-edge of the center mirror segment: 

Hres (1 — k(~ 1 j 2 1 ) \j ~ 1 %c j a (j represents the j th mirror segment counted from the 

center), at which point the mirror strength y j= '=0.7 increases with j. Fig.|2a) shows 7— j plot 
for different tapering profiles, which is the dependence of the filling fraction (/) on the index 
of mirror segment (J) starting from the center (j = 1). It can be seen that quadratically tapering 
profile results in linearly increasing mirror strengths, needed for Gaussian field attenuationQ]]. 
To verify this, we numerically solved the band diagram (Fig. Etb)) and obtained 7— / relation 
in Fig.|2c). As shown in Fig.Htd), linearly increasing mirror strength is indeed achieved after 
quadratic tapering. 



Next, with the optimized tapering profile, the cavity is formed by putting two such mirrors 
back to back, leaving a cavity length L in between (Fig. Q2 a ))- Fig- Be) shows the simulated 
g-factors for various Ls. Highest Qrad is achieved at zero cavity length (L=0), which supports 
the prediction in 0J based on ID model. 

Third, we verify that the cavity mode has a Gaussian-type attenuation profile. Fig. |2jg) 
shows the H z -fit\d distribution in the plane right above the cavity, obtained from 3D FDTD 
simulation. As shown in Fig. Efh), this field distribution can be ideally fitted with H z = 
sin(7Dc/a)exp(— <xt 2 )exp(— ^y 2 ), with a = 0.33, <7 = 0.14 and § = 14. The fitted value a agrees 
with the "period", and a agrees with that extracted value from Fig. |2d): a = = 0.13. 
Fig- Hi) shows H z distribution along the dashed line in Fig. [2jg)&(h). Therefore, we conclude 
that zero cavity length, fixed periodicity and a quadratic tapering of the filling fraction results 
in a Gaussian field profile, which leads to a high-Q cavity U. 

Finally, as we have pointed out in (T|, current method results in a cavity whose resonance is 
asymptotically approaching the dielectric band-edge frequency of the central mirror segment 
(circled in Fig. [2jb)). The deviation from the band-edge frequency can be calculated using 
perturbation theory EH El : 

SX_ = Jj#-j(rjDj^ 

A 2je\E\ 2 dV ( } 

E|| is the component of E that is parallel to the side wall surfaces of the holes and Dx is the 
component of D that is perpendicular to the side wall surfaces of the holes. Under Gaussian 
distribution, the major field component D y = cos(7r/ax)exp(— ax 2 )exp(— <i;y 2 ), 8s perturba- 
tion occurs at x = ±(j + 1 /2)a±Rj, where Rj = \J fjab/n denotes the radius of the f" hole 
(counted the center), with j=\,2...N, N is the total number of mirror segments at each side. 
Since the cavity mode has a Gaussian profile, 1 / \fa characterizes the effective length of the 
cavity mode, and scales linearly with N, with a nonzero intercept due to diffraction limit. For 
large N, the intercept can be neglected, and thus = \/2Q x 0.14 2 /N. Plug the perturbation 
induced by the quadratic tapering from / = 0.2 to / = 0.1 into Eqn. [6] the frequency offset 
5X /X v.sN can be obtained. Fig. |2jf) shows the frequency offset for different total number of 
mirror pairs (AO, calculated from the perturbation theory, as well as using FDTD simulations. It 
can be seen that the deviation decreases as the number of modulated mirror segments increases, 
and is below 1% for N > 15. 

Therefore, we verified that an ultrahigh-g, dielectric-mode cavity resonant at a target fre- 
quency can be designed using the following algorithm: 

(i) Determine a target frequency. For example in our case we want /target = 200THz. Since the 
cavity resonance is typically 1 % smaller than the dielectric band-edge of the central segment, 
estimated using the perturbation theory, we shift-up the target frequency by 1%, i.e. /adjusted = 
202THz. 

(ii) Pick the thickness of the nanobeam: this is often pre-determined by the choice of the 
wafer. For example, in our case, the thickness of the nanobeam is 220nm, determined by the 
thickness of the device layer of our silicon-on-insulator (SOI) wafer. 

(iii) Choose periodicity according to a = Xo/2n e g, where n e g is effective mode index of 
the cavity and can be estimated by numerical modeling of a strip waveguide that nanobeam 
cavity is based on. However, we found that the absolute value of the periodicity is not crucial in 
our design, as long as there exists a bandgap. Therefore, we pick « e ff = 2.23, which is a median 
value of possible effective indices in the case of free standing silicon nanobeam (n e g E [1 , 3.46]). 
This results in a = 330nm. 

(iv) Set nanobeam width. Large width increases the effective index of the cavity mode, pulls 
the mode away from the light line, thus reducing the in plane radiation loss. On the other hand, a 
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Fig. 3. (a) Schematic of the Gaussian nanobeam cavity, with an air hole in the symmetry 
plane (dashed line), (b) Energy distribution in the middle plane of the cavity obtained from 
3D FDTD simulation. (c)&(d) H z field distribution on the surface right above the cavity: 
(c) is obtained from 3D FDTD simulation and (d) is obtained from the analytical formula 
H : =cos(fx)exp(-CTx 2 )exp(-^j 2 ), with a = 0.33jUm,<T = 0.14, § = 14. (e) H z field dis- 
tribution along the dashed line in (c)&(d). Length unit: fim. 



large beam width will allow for higher order modes with the same symmetry as the fundamental 
mode of interest. Using band diagram simulations, we found that the width of 700nm is good 
trade-off between these two conditions (Fig. [2jb)). 

(v) Set the filling fraction of the first mirror section such that its dielectric band-edge is at the 
adjusted frequency: 202THz in our case. Band diagram calculations based on single unit cells 
are sufficient for this analysis. We found that an optimal filling fraction in our case is / start = 0.2 

(Fig. Kb)). 

(vi) Find the filling fraction that produces the maximum mirror strength for the target fre- 
quency. This involves calculating the mirror strength for several filling fractions (Fig.[2|c)), each 
of which takes one or two minutes on a laptop computer. In our case we found that / en( j = 0.1. 

(vii) Pick the number of mirror segments (AO to construct the Gaussian mirror: we found that 
N > 15 (on each side) are generally good to achieve high radiation-gs. 

(viii) Create the Gaussian mirror by tapering the filling fractions quadratically from /start 
(=0.2 in our case) to / en d (=0.1) over the period of N segments. From the above analysis, the 
mirror strengths can be linearized through quadric tapering (Fig.[2|d)). 

(ix) Finally, the cavity is formed by putting two Gaussian mirrors back to back, with no 
additional cavity length in between (L = 0). To achieve a radiation-limited cavity (2wg >> 
2radX 10 additional mirrors with the maximum mirror strength are placed on both ends of the 
Gaussian mirror. We will show in the next section, no additional mirrors are needed to achieve 
a waveguide-coupled cavity (2rad >> Gwg)- 

Besides the structure that were proposed in Q] (Fig. Q]), the alternative structure which has 
the air-hole in the symmetry plane, as shown in Fig. [3 a), also satisfies (i)-(ix). Both structures 
result in dielectric-mode cavities, since the bandgap of each mirror segments red-shifts away 
from the center of the cavity, and thus a potential well is created for the dielectric band-edge 
mode of the central segment. The difference is that the energy maximum in the air-hole centered 



cavity is no longer located in the middle of the structure, but instead in the dielectric region next 
to the central hole (Fig. Ob)). Fig.[3]c) shows the H z field profile in the plane right above the 
cavity, obtained from FDTD simulation. Fig. [3d) shows the fitted field profile using the exactly 
the same parameters that were used in the original structure shown in Fig.|2g)-(i), but with sine 
function replaced by cosine function. Fig.[3te) shows the H- distribution along the dashed line 
inFig.Hc)&(d). 

Armed with the analytical field profile of the cavities: H" dd (x) = sin(7Dc/a)exp(— Ox 2 ) and 
m ve "(x) = cos(nx/a)exp(-ax 2 ), we can obtain the radiation losses and and far fields of the 
cavities using the Fourier space analysis (33] ■ The Fourier transforms can be analytically ob- 
tained FT(H° dd ) = (exp(-(fc + 7r/a) 2 /4(j)-exp(-(fc-7r/fl) 2 /4a))//\/8o : andFT(//! :ven ) = 
(exp(— (£+ n/a) 2 /4a) - exp(-(fc - n / a) 2 / Aa)) / \/%a . Under aa 2 « 1, both distributions 
have their Fourier components strongly localized at k = ±n/a, as is verified by FDTD simula- 
tions in Fig.|4ja&b). Since H? (x) is an odd function, it always has a zero Fourier component 
at k = 0. Therefore, dielectric-centered cavities (Fig. [T} should have higher Q-factors. How- 
ever, in high-Q cavity designs, aa 2 << 1 is satisfied and thus both dielectric-centered and 
air-centered cavities have comparable Q-factors. FDTD simulation shows that the above H° dd 
and Hf m cavities have Qtot = 3.8 x 10 s and Qtot = 3.5 x 10 8 respectively. The mode volume 
of the ~H° dd cavity is 0.67(A res /n S i) 3 , smaller than the Hf en cavity (V = 0.76(Ar es /«si) 3 )- 

The far field radiation patterns (obtained using FDTD simulations) of the two cavities are 
shown in Fig.[4]c)&(d). The powers, in both cases, are radiated at shallow angles (> 70° zenith 
angle) to the direction of the waveguide. The H odd cavity has even less radiated power at small 
zenith angles, consistent with the above analysis. By integrating the zenith and azimuth angle 
dependent far field emission, we found that that 32% and 63% of the power emitted to +z 
direction can be collected by a NA=0.95 lens, respectively for H" dd cavity and H even cavity. 

3. Ultra-high Q, Dielectric-mode Photonic Crystal Nanobeam Cavities 

3.1. Radiation-Q limited and waveguide-coupled cavities 

Since the dielectric-centered H° cavity has smaller V than the Hf' en one, we focus our dis- 
cussion in the H odd case. Using the above design algorithm, we designed the Gaussian mirror 
and put 10 additional mirrors with the maximum mirror strength on both ends of the Gaussian 
mirror to obtain the radiation-limited cavity (Q W g >> Gtot)- We find in Fig.[5]that gtot increases 
exponentially and V increases linearly to the total number of mirror pairs in the Gaussian mirror 
(AO. A record ultra-high Q of 5.0 x 10 9 was achieved while maintaining the small mode volume 
of 0.9 x (Ar es /«si) 3 at AT = 30. 

Our design strategy has an additional important advantage over other types of photonic crys- 
tal cavities |[T7l - ll30l . that is: the cavity naturally couples to the feeding waveguide, as the hole 
radii decrease away from the center of the cavity. High-Q and high transmissions (T) cavities 
are possible with the above design steps (i)-(ix), with no additional "coupling sections" needed. 
We study T and dotal dependence on the total number of mirror pair segments in the Gaussian 
mirror (AO in Fig.|5jb). Partial g-factors (Qrad,2wg) were obtained from FDTD simulations, 
and T was obtained using T = 2totai/GwgEU- As shown in Fig.|5jb), we achieved a nanobeam 
cavity with Q = 1.3 x 10 7 ,T = 97% atN = 25. 

3.2. Higher order modes of the dielectric-mode cavity 

The ultra-high Q mode that we deterministically designed is the fundamental mode of the 
cavity. Meanwhile, higher order cavity modes also exist. The number of higher order modes 
depends on the width of the photonic band gap and total number of mirror segments in the 
Gaussian mirror. Fig. |6ja) shows the transmission spectrum of a waveguide-coupled nanobeam 
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Fig. 4. (a)&(b) The distribution of the spatial Fourier components of the cavity mode, ob- 
tained from 3D FDTD simulation: (a) for the H° cavity and (b) for the Ht ve " cavity 
respectively. (c)&(d) The far field profile of the cavity mode obtained from 3D FDTD sim- 
ulation: (c) for the H'J dd cavity and (d) for the Hf en cavity respectively. The inset cavity 
structure shows the orientation of the waveguide direction in (c)&(d). Dashed line indicates 
the symmetry plane. 
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Fig. 5. (a) Total g-factors (log(10) scale) and effective mode volumes (V/(Ar es /«si) 3 ) °f 
nanobeam cavities for different total number of mirror pair segments in the Gaussian mirror. 
In each case, 10 additional mirror segments with f=0.1 (maximum mirror strength) are 
added on both ends of the Gaussian mirror. Therefore, the total-g of the cavity is limited 
by radiation-g. A record ultra-high g of 5.0 x 10 9 is achieved with a Gaussian mirror that 
comprises 30 mirror segments and an additional 10 mirror pairs on both ends, (b) On- 
resonance transmissions and total g-factors (log(10) scale) v.s the total number of mirror 
pair segments in the Gaussian mirror. In this case additional mirror pairs (10 of them) are 
not included. A record high-T (97%) and high-g (1.3 x 10 7 ) cavity is achieved at N = 25. 

cavity, that has 12 mirror pair segments in the Gaussian mirror. Since the cavity was designed 
to be a waveguide-coupled cavity, the simulation was performed by exciting the input waveg- 
uide with a waveguide mode, and monitoring the transmission through the cavity at the output 
waveguide. The band-edge modes are observed at wavelengths longer than 1.6/im and shorter 
than 1.3jtim. Fig.[6jb)-(d) shows the major field-component (E y ) distribution of the three cavity 
modes. As expected, the eigenmodes alternate between symmetric and anti-symmetric modes. 
Symmetry plane is defined perpendicular to the beam direction, in the middle of the cavity 
(dashed line in Fig. [(J. We note that transversely odd modes are well separated from the trans- 
versely symmetric cavity modes, hence were not considered in Fig. [6] 

4. Ultra-high Q, Air-mode Photonic Crystal Nanobeam Cavities 

4.1. Radiation-Q limited cavity 

An air-mode cavity concentrates the optical energy in the low index region of the cavity. There- 
fore, these cavities are of interest for applications where strong interactions between light 
and material placed in the low index region of the cavity is required, including nonlinear 
enhancement^ |, optical trappingQ, biochemical sensing|8| and light-atom interactional]. 
The ultra-high Q air-mode nanobeam cavity is realized by pulling the air-band mode of pho- 
tonic crystal into its bandgap, which can also be designed using the same design principles that 
we developed for dielectric-mode cavities. In contrast to the dielectric-mode case, the resonant 
frequency of the air-mode cavity is determined by the air band-edge frequency of the center 
mirror segment. And to create the Gaussian confinement, the bandgaps of the mirror segments 
should shift to higher frequencies as they go away from the center of the cavity. This can be 
achieved by progressively increasing the filling fractions of the mirror segments away from the 
center of the structure (instead of decreasing in the dielectric -mode cavity case). One way to 
accomplish this is to increase the size of the holes away from the center of the cavity. While 
this may be suitable for non-waveguide coupled (radiation-Q limited) cavities, it is not ideal for 
a waveguide-coupled cavity, where high transmission efficiency through the cavity is required. 
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Fig. 6. (a) Transmission spectrum of the cavity from FDTD simulation, (b)-(d) The E y 
field distribution in the middle plain of the nanobeam cavity. Resonances and symmetries 
of the modes are indicated in the plot. Symmetry plane is perpendicular to the nanobeam 
direction, in the middle of the cavity, as indicated by the dashed line. Length unit in (b)-(d) 
is ^m. 




Fig. 7. (a) TE band diagram for an air-mode nanobeam cavity. Hole radii r = lOOnm.a = 
330nm, b=l/im (red) and b=0.7/lm (black), (b) Mirror strengths for different beam widths, 
(c) Linearization of mirror strengths after quadratic tapering the beam widths. 
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Fig. 8. (a) Total g-factors (log(10) scale) and effective mode volumes (V/(Ar es /«si) 3 ) °f 
the nanobeam cavities for different total number of mirror pair segments in the Gaussian 
mirror. In each case, 10 additional mirror segments with w=0.7/lm are added on both ends 
of the Gaussian mirror, so that the total-g of the cavity is limited by radiation-g. A record 
ultra-high g of 1.4 x 10 9 is achieved with a Gaussian mirror that comprises 30 mirror 
segments and 10 additional mirror pairs on both ends, (b) On-resonance transmissions and 
total g-factors (log(10) scale) v.s the total number of mirror pair segments in the Gaussian 
mirror. In this case additional mirror pairs (10 of them) are not included. A record high-J 
(96%) and high-g (3.0 x 10 6 ) cavity is achieved at N = 25. 

For this reason, we employ the design that relies on tapering of the waveguide width instead 
of the hole size. Similar geometry was recently proposed by Ahn et. al. l42l for the design of a 
dielectric-mode photonic crystal laser with different design method. 

The same design steps can be followed as in the dielectric-mode cavity case, with the fol- 
lowing changes: First, the adjusted frequency (198THz) is 1% lower than the target frequency 
(200THz). (The thickness of the nanobeam is 220nm and period is 330nm, same as previous 
case.) Second, the nanobeam width at the center of the cavity is w s tart = 1/im (Fig. |7ja)), with 
the hole radii kept constant at lOOnm. Third, to create the Gaussian mirror, the beam widths 
are quadratically tapered from w start = Ifim to w en d = 0.7/im, which produces the maximum 
mirror strength (band diagrams shown in Fig. [jj a)). This procedure involves calculating the 
mirror strength for several beam widths (Fig. [71b)), each takes one or two minutes on a laptop 
computer. As shown in Fig. |7jc), the mirror strengths are linearized after the quadratic taper- 
ing. In order to achieve a radiation-g limited cavity, 10 additional mirror segments are placed 
at both ends of the Gaussian mirror that has beam width w cn( j = 0.7/im. 

Similar in the dielectric-mode cavity cases, H° dd and H e , ven air-mode cavities can be formed 
by placing the air and dielectric in the central symmetric plane of the cavity, respectively. Again, 
we will focus on H° dd , air-mode cavities and the conclusions will be valid to the H- ve " cavities 
as well. Fig. [Ua) shows the total Q of nanobeam cavities hat have different total number of 
mirror pair segments in the Gaussian mirrors. We have achieved a record ultra-high Q of 1 .4 x 
10 9 , air-mode nanobeam cavity. As shown in Fig. [HI a), the effective mode volumes of the air- 
mode cavities are much larger than the dielectric-mode cavities. 

4.2. Cavity strongly coupled to the feeding waveguide 

As we have pointed out, the tapering-width approach (as compared to taping hole radii) offers 
a natural way of coupling the nanobeam air-mode cavity to the feeding waveguide. Since the 
width of the beam is decreasing, the cavity naturally couples to the feeding waveguide. We 
study T and gtotai dependence on the total number of mirror pair segments in the Gaussian 
mirror (AO using FDTD simulations. As shown in Fig.[8jb), we were able to design a nanobeam 
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Fig. 9. (a) Transmission spectrum of the cavity from FDTD simulation. (b)&(c) The E y 
field distribution in the middle plain of the nanobeam cavity. Resonances and symmetries 
of the modes are indicated in the plot. Symmetry plane is perpendicular to the nanobeam 
direction, in the middle of the cavity, as indicated by the dashed line. Length unit in (b)&(c) 
is flm. 



cavity with Q = 3.0 x 10 6 ,T = 96% at N = 25. 
4.3. Higher order modes of the air-mode cavity 

The ultra-high Q cavity that we were able to design is the fundamental mode of the cavity. 
Higher order modes coexist with the fundamental modes inside the band gap. Fig. |9ja) shows 
the transmission spectrum of a waveguide-coupled air-mode nanobeam cavity, that has 15 mir- 
ror pair segments in the Gaussian mirror. The band-edge modes are observed at wavelengths 
longer than 1 .6/im. The modes in the range of 1 .2/im to 1 .35/xm are formed by the higher order 
band modes in Fig. |7ja). Fig.|9jb)-(c) shows the major field-component distribution (Ey) of the 
two cavity modes inside the bandgap. 

5. Conclusion 

We have presented a detailed analysis and a deterministic design of the ultra-high Q pho- 
tonic crystal nanobeam cavities. With this method, Q > 10 9 radiation-limited cavity, and 
Q > 10 7 , T > 95% waveguide-coupled cavity were deterministically designed. These g-factors 
are comparable with those found in whispering gallery mode (WGM) cavities l43l - ll45l and 
have ultra-small mode volumes, typically two or three orders of magnitude smaller than WGM 
ones. Furthermore, energy maximum can be localized in either the dielectric region or air re- 
gion with this method. Although we demonstrated designs for TE-like, transversely symmetric 
cavity modes, the design method is universal, and can be applied to realize nanobeam cavi- 
ties that support TM-polarized modes, as well as line-defect 2D photonic crystal cavities. We 
believe that the proposed method will greatly ease the processes of high Q nanobeam cavity 



design, and thus enable both fundamental studies in strong light and matter interactions, and 
practical applications in novel light sources, functional optical components (filters, delay lines, 
sensors) and densely integrated photonic circuits. 

6. Acknowledgments 



We acknowledge numerous fruitful discussions with M.W. MuCutcheon and P. B. Deotare. This 
work is supported by NSF Grant No. ECCS-0701417 and NSF CAREER grant. 



